Quantum Nature of Relaxation of Paramagnetic and Optical Systems by Strong 
Dipole-Photon and Dipole-Phonon Coupling 
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Matrix-operator difference-differential equations for dynamics of spectroscopic transitions in ID 
multiqubit exchange coupled (para)magnetic and optical systems by strong dipole-photon and 
dipole-plionon coupling are derived within the framework of quantum electrodynamics and quantum 
phonon field theory. It has been established, that in the model considered the relaxation processes 
are of pure quantum character, which is determined by the formation of the coherent system of 
the resonance phonons and by the appearence along with absorption process of EM-field energy 
the coherent emission process, acompanying by phonon Rabi quantum oscillation, which can be 
time-shared. For the case of radiospectroscopy it corresponds to the possibility of the simultaneous 
observation along with (para)magntic spin resonance the acoustic spin resonance. Theoretical con- 
clusions are confirmed experimentally in radiospectroscopy. It has been found in particular, that the 
lifetime of coherent state of collective subsystem of resonance phonons in disordered carbon samples 
- anthracites of medium-scale metamorphism - is very long and even by room temperature it is 
evaluated in ~ 10 ms. The phenomenon of the formation of the coherent system of the resonance 
phonons can be used in a number of practical applications, in particular by elaboration of logic 
quantum systems including quantum computers and quantum communication systems. 

PACS numbers: 42.50.Ct, 61.46.Fg, 73.22.-f, 78.67.Ch, 77.90.+k, 76.50.-l-g 



I. INTRODUCTION AND BACKGROUND 

Quantum electrodynamics (QED) and quantum field 
theory at all take on more and more significance for its 
practical application and thy, in fact, become to be work- 
ing instrument in spectroscopy studies and industrial 
spectroscopy control. Quantum dynamics of two level 
systems (qubits), coupled to a single mode of an electro- 
magnetic cavity are of considerable interest in connection 
with quickly devefoping new quantum physics branches 
like to cavity quantum electrodynamics [l| or quantum 
computing Spin-based solid state quantum bits are 
known to have long coherence times, while also offering 
the promise of scalability, and are natural building blocks 
for quantum computation. Phosphorus donor nuclei in 
silicon have been known since the 1950s to have some 
of the best spin coherence properties in solids. The spin 
coherence time T2 measured by Hahn spin echo method 
for donor electron spins in bulk Si:P has been reported 
to be equal ~ 60ms [1]. This is the longest coherence 
time measured in electron spin qubits, and greatly ex- 
ceeds the value reported, for instance, in GaAs quantum 
dots, which, measured also by Hahn spin echo method, 
is equal to ~ 100/j,s I^]. However, fabrication of ordered 
and gated donor arrays and coherent control over donor 
electrons has turned out to be extremely difficult. There 
is now much interest in fabricating spin-based devices in 
diamond, with potential applications in quantum com- 
munication, quantum computation, and magnetometry. 
Nitrogen-vacancy (NV) centers located deep within a di- 
amond lattice appear promising solid-state spin qubits 
since they combine optical initialization and readout ca- 



pabilities owing to long electron spin coherence life times 
approaching 1ms at room temperature Q , and the abil- 
ity to control coupling to individual nuclear spins. In all 
of aforesaid applications it is necessary or at least advan- 
tageous to couple NV centers to optical structures like to 
waveguides and microresonators, to enable communica- 
tion between distant qubits or to allow efficient extrac- 
tion of emitted photons. Therefore, a reliable method is 
needed to create NV centers with good spectral proper- 
ties in close proximity (< lOOnm) to a diamond surface. 
In addition, the charge state of the NV center must be 
controlled. Given tasks are under study at present. We 
will show in present work another way to obtain long- 
lived coherent states with similar field of practical ap- 
plication. The prediction will be based on QED-theory 
of spectroscopic transitions. The simplest models which 
capture the salient features of the relevant physics in 
given field are the Jaynes-Cummings model (JCM) 6] for 
the one qubit case and its generalization for multiqubit 
systems by Tavis and Cummings Tavis-Cummings 
model was generalized in by taking into account the 
ID-coupling between qubits. Recently QED-model for 
one chain coupled qubit system was generalized for qua- 
sionedimensional axially symmetric multichain coupled 
qubit system It is substantial, that in the model, 
proposed in Q the interaction of quantized EM-field with 
multichain qubit system is considered by taking into ac- 
count both the intrachain and interchain qubit coupling 
without restriction on their values. It follows from theo- 
retical results in , P and from their experimental con- 
firmation in [1^, that by strong interaction of EM-field 
with matter the correct description of spectroscopic tran- 
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sitions including stationary spectroscopy is achieved the 
only in the frame of QED consideration. It concerns both 
optical and radio spectroscopies, that means, that QED 
consideration has to be also undertaken by electron spin 
resonance (ESR) studies in the case of strong interaction 
of EM-field with spin systems. It is reasonable to sug- 
gest, that analogous conclusion can be drawn for the case 
of strong interaction of phonon field with spin system or 
electron system. In other words it seems to be reasonable 
the idea, that relaxation of paramagnetic (or optical) cen- 
ters in the case of strong spin-phonon (electron-phonon) 
interaction can be described correctly the only in the 
frames of quantum field theory. 

The aim of the work presented is to derive the system of 
equations for dynamics of spectroscopic transitions in ID 
multiqubit exchange coupled (para) magnetic and optical 
systems by strong dipole-photon and dipole-phonon cou- 
pling within the framework of quantum electrodynamics 
and quantum phonon field theory (QPTT) and to show 
both theoretically and experimentally, that new quantum 
plysics phenomenon - the formation of longlived coher- 
ent state of resonance phonons leading to appearance uf 
quantum acoustic Rabi oscillations takes place and to 
propose the area of its application. 
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II. RESULTS AND DISCUSSION 

A. Quantum Equations for Spectroscopic 
Transitions by Strong Electron-Photon and 
Electron-Phonon Coupling 

Recently in the work the system of difference- 
differencial equations for dynamics of spectroscopic tran- 
sitions for both radio- and optical spectroscopy for the 
model, representing itself the ID-chain of N two-level 
equivalent elements coupled by exchange interaction (or 
its optical analogue for the optical transitions) betweenn 
themselves and interacting with quantized EM-field and 
quantized phonon field its optical has recently been de- 
rived. Naturally the equations are true for any 3D sys- 
tem of paramagnetic centers (PC) or optical centers by 
the absence of exchange interaction. In given case the 
model presented differs from Tavis-Cummings model 
by inclusion into consideration of quantized phonon sys- 
tem, describing the relaxation processes from quantum 
fied theory position. Seven equations for the seven oper- 
ator variables, describing joint system {field -I- matter} 
can be presented in matrix form by three matrix equa- 
tions. They are the following 
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is vector-operator of spectroscopic transitions for /th 
chain unit, I = 2,iV — 1 Its components, that is, 

the operators 
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are set up in correspondence to the states |jt,),(TO„|, 
where v — 1,N, j = a^(3^ m — a, 13. For instance, the 
relationships for commutation rules are 
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where vector operators I — 2, N — 1, are given 

by the expressions 

Gi-i,i+i = Gili^i+ie+ + i_,+ie_ + ^f_i_j+ie;, (8) 
in which 

k 







'Fr' 










= 2||g|| 




(1) 


k 


- ^! - 




.Ft . 






Gi-i,i+i = -^i - 



I' 



'l+l 



(9b) 



(9c) 



3 



Here operator f^j: is 

fit: 
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In relations ([S]) J is the exchange interaction constant in 
the case of magnetic resonance transitions or its optical 
analogue in the case of optical transitions, the function 
Vik in (0 is 
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where pj™ is matrix element of operator of magnetic 
(electric) dipole moment Pi of l-th chain unit between 
the states and \mi) with j e {q;,/3}, m € {a, 13}, 
j 7^ TO, is unit polarization vector, ep^ is unit vec- 
tor along P/-direction, £^ is the quantity, which has the 

dimension of magnetic (electric) field strength, k is quan- 
tized EM-field wave vector, the components of which get 
a discrete set of values, ujj: is the frequency, corresponding 
to fcth mode of EM-field, ai and dr are EM-field creation 

k ^ 

and annihilation operators correspondingly. In the sug- 
gestion, that the contribution of spontaneous emission is 



relatively small, we will have p[ 
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Here B^^^ is 
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hi (hq) is the creation (annihilation) operator of the 
phonon with impulse q and with energy hw^, is 
electron-phonon coupling constant. In equations ^ and 
([3]) llfTpll is Pauli z-matrix, in equation ([!]) is diagonal 
matrix, numerical values of its elements are dependent on 
the basis choice. It is at appropriate basis 
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Right hand side expression in ([T]) is vector product 
of vector operators. It can be calculated by using of 
known expression ((T5|) with additional coefficient 
only, which is appeared, since 
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the products of two components of two vector operators 
are replaced by anticommutators of corresponding com- 
ponents. Given detail is mapped by symbol (g) in ([71) and 
by symbol ' in determinant (jlSp . 

It follows from comparison with semiclassical Landau- 
Lifshitz (L-L) equation for dynamics of spectroscopic 
transitions in a chain of exchange coupled centers [l5|, 
p^ . that the equation, which is given by ([T]) is its 
QED-generalization. In comparison with semiclassical 
description, where the description of dynamics of spec- 
troscopic transitions is exhausted by one vector equa- 
tion (L-L equation or L-L based equation), in the case 
of completely quantum consideration L-L type equation 
describes the only one subsystem of three-part-system, 
which consist of EM-field, dipole moments' (magnetic or 
electric) matter subsystem and phonon subsystem. It 
was concluded in [1J|, that the presence of additional 
equations for description of transition dynamics by QED 
model in comparison with semiclassical model leads to a 
number of new effects, which can be predicted the only 
by QED consideration of resonance transition phenom- 
ena. One of new effect was described in [l3|, starting 
the only from the mathematical structure of the equa- 
tions. It was argued, that the equations ([T]), ^ represent 
themselves vector-operator diffcrcncc-difi'erential gener- 
alization of the system, which belongs to well known fam- 
ily of equation systems - Volterra model systems, widely 
used in biological tasks of population dynamics studies, 
which in its turn is generalization of Verhulst equation. 
In other words, it was predicted, for instance, that by 
some parameters in two-sybsystem Volterra model the 
stochastic component in solution will be appeared. Given 
prediction has experimental confirmation by the study 
of optical properties in carbynes fl6j indicating, that in 
given material strong electron-photon interaction is re- 
alized, which allows to explain the posibility to observe 
the stationary IR-reflection or absorbtion spectra both in 
usual and in stochastic regime. 

The terms like to right hand side terms in (jS]) were 
used in so called "spin-boson" Hamiltonian ,17] and in 
so called " independent boson model" [ISj] . Given models 
were used to study phonon effects in a single quantum 
dot within a microcavity fl9ll . (20j . [2l| . [23 |. [23j . So, 
it has been shown in 1231 , [2j] , that the presence of the 
term in Hamiltonian 1141 
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which coincides with corresponding term in Hamiltonian 
in [l^], [l^l at = 1 [contribution of given term to 

N 

the equations for spectroscopic transitions is ±^ ^i^qy 

1=1 

see equation (note that the equations for spec- 

troscopic transitions were not derived in above cited 
works [3, d^, [U, [13, HI)] leads the only to ex- 
ponential decrease of the magnitude of quantum Rabi 
oscillations with increase of electron-phonon coupling 
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strength and even to their supression at relatively strong 
electron-phonon coupling. However by strong electron- 
photon coupling amd strong electron-phonon coupling 
quite other picture of quantum relaxation processes be- 
comes to be possible. Really, if to define the wave func- 
tion of the chain system, interactig with quantized EM- 
field and with quantized lattice vibration field, to be vec- 
tor of the state in Hilbert space over quaternion ring, that 
is quat ernion function of quaternion argument, then like 
to [ij] can be shown, that the equations (P) to ([3]) are 
Lorentz invariant and the transfer to observables can be 
realized. In particular, taking into account, that quater- 
nion vector of the state is proportional to spin, the Hamil- 
tonian, given by (|16|) describes in fact the interaction of 
phonon field with z- component of the spin of matter 
subsystem. It seems to be reasonable to take into consid- 
eration the interaction of phonon field with S^- and 
components of the spin of matter subsystem. Therefore 
we come in a natural way to the following Hamiltonian 
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where TiP \s chain Hamiltonian by the absence of the in- 
teraction with EM-field, "H^ is field Hamiltonian, ikP^ 
is Hamiltonian, describing the interaction between quan- 
tized EM-field and atomic chain. Hamiltonian is 



(18) 



where w' is chain Hamiltonian in the absence of the inter- 
action between structural elementary units of the chain. 
■H" is given by the expression 
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Here m — a,P, Emv are eigenvalues of 'W' , which corre- 
spond to the states \m^) of vth chain unit. Hamiltonian 
is 

(20) 

N ^ ' 



It is suggested in the model, that and |/3„) are eigen- 
states, producing the full set for each of N elements. It 
is evident, that given assumption can be realized strictly 
the only by the absence of the interaction between the 
elements. At the same time proposed model will rather 
well describe the real case, if the interaction energy of ad- 
jacent elements is much less of the energy of the splitting 
^0 = — £a between the energy levels, corresponding 
to the states |a„) and |/3„). The case considered includes 
in fact all known experimental situations. It is clear, 
that Hamiltonian H'-^^ of interaction of quantized EM- 
field with atomic chains can also be represented in the set 
of variables, which includes the components of spectro- 
scopic transition vector operator ay. Really, suggesting 



dipole approximation to be true and polarization of field 
components to be fixed, we have 
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where p^™ is matrix element of operator of magnetic 

(electric) dipole moment Pj of j-th chain unit between 
the states and \mj) with Ij — aj,Pj,mj = aj,f3j, 
is unit polarization vector, ep is unit vector along Pj- 
dircction, 2^: is the quantity, which has the dimension of 

magnetic (electric) field strength, k is wave vector, a-j^ is 
field annihilation operator. In the suggestion, that the 
contribution of spontaneous emission is relatively small, 
we will have p^™ = p™' = pj, where I — a, /3,m — a, (3. 
Let define the function 
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Then the expression (1211) can be rewritten in the form 

^"'^ -jZTMjki^J+-t)H+^J+-TH^d^ (23) 



where at is EM-field creation operator, a- is EM-field 

k k 

annihilation operator, superscript * in q* jj: means com- 
plex conjugation. Field Hamiltonians have usual form 



for EM-field and 



for phonon field. Hamiltonian V. is 

n:,CPh _ ^CPh I -£,CPh 
it — rt^ -r ti-j^ , 

where H^^'^ is determined by the expression 
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Hamiltonian H^^^ can be represented in the following 
form 

= E E ^^(-7 + (Ai)*(^- + ^m- 

(28) 



and A± ^ 

which characterisire correspondingly the interaction with 



Here and A? are electron-phonon coupling constants, 
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z- component 5*1 and with 5+- and S~ components of 
the spin of jth chain unit. It seems to be understand- 
able, that they can be different in general case. More- 
over, in order to take into account the interaction with 
both equilibrium and noncquilibrium phonons both the 
electron-phonon coupling constants have to be complex 
numbers, that takes proper account by expressions (1271) . 



It can be shown, that the equations of the motion for 
spectroscopic transition operators (?;, for quantized EM- 
field operators a^, at and for phonon field operators 

&i are the following. Instead equation ([Ij we have 
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The equation ^ remains without changes. The equation 
® is 
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Thus, QED - QPFT model for dynamics of spectro- 
scopic transitions in ID multiqubit exchange coupled sys- 
tem is generalized by taking into account the earlier proof 
[13 ] , that spin vector is quaternion vector of the state of 
any quantum systen in Hilbert space defined over quater- 
nion ring and consequently all the spin components has 
to be taken into account. New quantum phenomenon 
is predicted. The prediction results from the structure 
of the equations derived and it consists in the following. 
The coherent system of the resonance phonons, that is, 
the phonons with the energy, egualed to resonance pho- 
ton energy can be formed by resonance, that can lead to 
appearance along with Rabi oscillations determined by 
spin (electron)-photon coupling with the frequency ft^^ 
of Rabi oscillations determined by spin (electron)-phonon 
coupling with the frequency fi^^''. In other words QED 
- QPFT model predicts the oscillation character of quan- 
tum relaxation, that is quite different character in com- 
parison with phenomenological and semiclassical Bloch 
models. Moreover if | \< g the second Rabi oscilla- 
tion process will be observed by stationary state of two 
subsystems { EM-fied -|- magnetic (electric) dipoles }, 
that is, it will be registered in quadrature with the first 
Rabi oscillation process. It can be experimentally de- 
tected even by stationary spectroscopy methods. 



B. Experiment 

ESR studies in anthracite samples of medium-scale 
metamorphism have been carried out on stationary ESR 
spectrometer "Radiopan" at room temperature. The 
100 kHz high frequency modulation of static magnetic 
field Hq was used. The static magnetic field Hq was 
sweeped in two directions from lesser value to bigger 
value (^^ > 0) and from bigger value to lesser value 
(^^ < 0). Two regime of registration were used - with 
automatic microwave frequency adjustment (AF-mode) 
and without automatic microwave frequency adjustment 
with constant microwave frequency (CF-mode). The 
spectra obtained are presented in Figures 1, 3-9. The 
most interesting, that the spectra can be registered by 
using of 100 kHz high frequency (HE) modulation canal 
at modulation amplitude Hm = 0, Figure 1. It is seen 
from Figure 1, that the spectrum has very unusual for 
stationary ESR-spectroscopy form, it is similar partly to 
the spectra, registered by nonstationary transient ESR- 
spectroscopy. At the same time the functional depen- 
dence is another. The spectral distribution of nutation 
signals in transient ESR-spectroscopy is described usu- 
ally by continuous Bessel functions. The spectrum, pre- 
sented in Figure 1, has time-discrete character and con- 
sist of two group of very narrow ESR lines. For instance, 
the linewidth AH of the central line with maximal am- 
plitude in left group is equal to 0.01 (±0.004) G. It allows 
to determine (/-value (that is g-factor in Zeeman term of 
spin-Hamiltonian) very precisely. It is equal to 2.002735 
(±0.000001). The linewidth of the other lines is prac- 
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Figure 1: Quantum Rabi oscillation process indicating on the 
formation of the coherent system of the resonance phonons 



tically the same. The lines presented in Figure 1 seem 
to be the most narrow among all the lines, registered 
in ESR-spectroscopy at all. The position of the central 
line with maximal amplitude in right group is given by g- 
value 2.00255 (±0.00001). The structure of the spectrum 
is very similar to theoretical time evolution for a system 
with one qubit in JCM, Figure 2. Really, the spectra in 
both the groups can be described by almost sinusoidal de- 
pendence with asymmetric Gaussian envelope. Figure 3. 
The qualitative difference between time evolution given 
by JCM and spectral distribution presented in Figure 1 
is that, that the left group of lines in JCM has cosinu- 
soidal dependence. Figure 2. The possibility to explain 
the spectrum observed by usual amplitude modulation 
from possible technical noise (determined for instance by 
the operation of static magnetic field stabilization unit or 
microwave frequency stabilization unit) is excluded, since 
the envelope by amplitude modulation is described by si- 
nusoidal dependence, we have the envelope, described by 
Gauss distribution with slightly different parameters for 
left and right sides. So the fit of experimental data. Fig- 
ure 3, for left hand side group of lines is given by the 
following expression 



y 2/0 + 
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(35) 



where yo = 52.4, He = 3322. 95G, W = 0.36G, A = 32G. 

Further the value of g-factor of the center of the first 
group agrees well with g-factor of ESR line rof the same 
sample registered by usual conditions with nonzero am- 
plitude of 100 kHz high frequency modulation [however 
it was determined substantionally more precisely]. The 
dependence of the amplitude and the shape of the first- 
group-lines and dependence of the amplitude, the shape 
and the position of the second-group-lines ( which are 
not analysed in given communication on the microwave 
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Figure 2: Theoretical evolution for JCM-system, initial state 
It/j) ® \a) of which is direct product of atomic ground state 
and coherent field state with n = 50 
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Figure 3: The fit of left hand side envelope of left oscillation 
group in the spectrum, presented in Figure 1 by Gaussian 



power) is addditional argument for reality of derection 
of spin system response at Hm = 0. Given arguments 
and comparison with Figure 2 allow to conclude that the 
spectrum in Figure 1 represents itself the quantum Rabi 
oscillation picture. Really, let us touch on JCM in some 
details. Central role in JCM plays the sum with infinite 
summation limit, which represents on a scale [-1,1] the 
degree of excitation of two-level system resulting of inter- 
action between the atom dipole or spin and single mode 
quantized EM-ficld. It is 

°° lap" 

(o-^(i)) = _exp[-^|ap] ^ J-^cos2gv^t, (36) 
n=i 

where \a) is fully coherent state of the field, taking place 
at i = 0, at that |ap = n, that is, it is the average 
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Figure 4: The fit of the splitting between the lines in left 
oscillation group in the spectrum, presented in Figure 1 by 
linear dependence 



Figure 6: ESR spectrun of anthracite sample by automatic 
microwave frequency adjustment, < 



6 

T3 

+^ 

"q. 

E 4 
< 




3351 3352 

Magnetic Field (G) 



3353 



"q. 
E 
< 




3352,26 3352,28 3352,30 3352,32 3352,34 
Magnetic Field (G) 



Figure 5: ESR spectrun of anthracite sample by automatic 
microwave frequency adjustment, > 



Figure 7: Central part of ESR spectrun of anthracite sample 



with automatic microwave frequency adjustment, 



> 



number of photons in the field, g is couphng constant 
between field and atom (spin). The sum in p6p cannot 
be expressed exactly in analytical form. For very short 
times and very large n behavior of ((7^(i)) is determined 
by cos2g^/^t. Cummings [l^] has shown, that by reso- 
nance and for intermediate time t values the cosine Rabi 
oscillation damp quickly (so called collapse takes place). 
Given damping can be described by Gaussian envelope 



exp[-^(5t)'] 



(37) 



It is substantial, that it not depends on field intensity 
unlike to semiclassical Rabi oscillation damping process 



and it is determined entirely the o nly by coupling con- 
stant g. The authors of the work 25| have found, that 
JCM contains so called revival process with revival time 
Tr, given by the expression 
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(38) 



where A is deviation of field mode frequency from res- 
onance value. Revival process takes place at all time 
values, satisfying the relation t = kTfj, k G N. It is 
seen from ([38|) that revival time depends on n and it is 
proportional to field amplitude at A = like to depen- 
dence of Rabi frequency by semiclassical consideration. 
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Figure 8: Central part of ESR spectrun of anthracite sample 
with automatic microwave frequency adjustment, < 



Let us illustrate the dynamics of spectroscopic transi- 
tions in JCM. The initial state for single qubit system 
can be defined by direct product of the two level matter 
(atomic, spin and so on) subsystem state jV') and coher- 
ent EM- field state \a). It is the following 



where IV') is 



and la) is 



l*(0)) = 



ClIV'l) + C2|V'2) 



exp[ 



n ,a 
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(39) 
(40) 

(41) 



The Rabi oscillations of the probability, that the qubit is 
in the initial state in the first damping stage (collapse), 
on a time scale of — — and then revive at Tr ~ 27rVn 
are illustrated in Figure 2 for ci = 1, C2 = by plotting 



^|(V'i,n|*(t))p 



(42) 



where {ipi,n\ = {ipil (8) (n| is the ground state for the 
qubit with n photons in the cavity. Our measurements 
of the dependence of the ESR spectral characteristics on 
the microwave power (photon number) is also showed 
that in correspondence with JCM the positions of the 
lines of left hand side group are not dependent on pho- 
ton number, at the same time the positions of the lines 
of right hand side group are dependent on photon num- 
ber. The aforesaid difference between our results and 
JCM can be easily explained if we determine the Rabi 
frequency value. It is equaled to 28 kHz. Given value is 
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substantially lesser in the same microwave power range, 
than the Rabi frequency values in excited by microwave 
power systems. So, for instance, in Si the Rabi frequency 



changes from sa 1 x 10°rads 



to 



4.5 X 10^rads~ 



for the Si-P3 centers [28|. It means that we really have 
observed the quantum relaxation oscillations, predicted 
theoretically in Section 2A. The prediction of sinusoidal 
character (instead cosinusoidal in JCM) is also confirmed. 
The comparison with the model for dynamic of spectro- 
scopic transitions, given in Section 2A is correct, since at 
modulation amplitude Hm — the absorption cannot be 
registered. At the same time the emission can be regis- 
tered, since by detection at 100 kHz the quantity which 
changes essentially slowly is registered like to usual ESR 
detection. Then, if to omit the part from Hamiltonian, 
corresponding to EM-field and z-relaxation term and to 
set N equaled to 1, we will have in rotating wave approx- 
imation JCM Hamiltonian for phonon subsystem. The 
linear dependence of splitting, presented in Figure 4, can 
be explained by the dependence of fl^^^ on the static 
H-field, since by the change of static H-field the distance 
between the energy levels in qubits is also changes. Ex- 
panding the dependence n^^^{H) in Taylor series and 
restricting by linear term (it is correct, since H-change is 
small), we obtain the agreement with experiment. 

To confirm the conclusion on coherent emission char- 
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acter (which is hke to maser emission character) of the 
spectrum, observed by = and to establish the struc- 
ture of the centers, which can emit stationary, we have 
undertaken the study of ESR-responce at various values 
of HF-modulation amplitude. The spectra, registered by 
100 kHz HF-modulation of static magnetic field with am- 
plitude Hm — O.OIG are presented in Figures 5 to 9. The 
lines have also unusual shape for absorption derivative. 
The spectra are also very different by the change of sweep 
direction, compare Figures 5 and 6, 7 and 8, which all 
were registered with automatic microwave frequency ad- 
justment. It is seen from Figures 5 and 6, that for both 
static field sweep directions one of the spectrum wings 
has stepwise shape, which usually is characteristic by AF 
adjustment cycle skip. At the same time cycle skip does 
not takes place, that is confirmed by detailed registra- 
tion of central part of the spectra in Figures 5 and 6, see 
Figures 7 and 8, and by simultaneous (parallel) registra- 
tion of dependence of the frequency of measuring cavity, 
that is time dispersion registration. Figure 9. It is seen 
from Figure 9, that microwave frequency jumps are not 
exceeding 1 MHz, that is, they are on the order of value 
less than the adjustment range of AF adjustment unit, 
which is ~ 10 MHz. It is seen also from Figure 9, that 
time dispersion responce has practically the same am- 
plitude for the both sweep directions in distinction from 
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the spetra, presented in Figures 5 and 6, where the am- 
phtude and intensity of resonance responce at > 
markendly exceed the amplitude and the intensity of res- 
onance responce at < 0. So, amphtude ratio is equal 
2, intensity ratio is evaluated to be equal ~ 4.6. It is in- 
teresting, that the step position is different for > 
and < 0, that is original hysteresis is taking place. 
All the differences for both sweep directions disappear by 
the registration without automatic microwave frequency 
adjustment, and for both sweep directions we obtain the 
spectrum, presented in Figure 10. At the same time the 
line has Dyson shape [l^ , which is characteristic for ESR 
absorbtion in metals and it is determined by the space 
dispersion contribution (27| . which is appeared in con- 
ductive media. Dispersionlike shape of absorption ESR- 
response is also the indication that the ESR centers are 
also strongly mobile [in the case of static PC in conduc- 
tive media the space dispersion contribution and absorp- 
tion contribution to resulting ESR responce is 1 to 1 [27]]. 
On Dyson effect with the strong mobility (strong spin dif- 
fusion) of PC indicates the evident deviation of line shape 
from Lorentz or Gauss usual dispersion line shapes. Re- 
ally, experimental value of asymmetry extent is equal 
to 23.4±0.5. (Compare with ■§ — & and ^ — 3.5 for usual 
Lorentz and Gauss dispersion line shapes correspond- 
ingly). It means in its turn, that in anthracite samples 
the electron-(spin)-photon interaction is really strong for 
both resonance and nonresonance interations. In result 
at relatively large amount of PC (necessary condition, 
which was formulated earlier) the sufficient condition for 
formation of the second cavity in measuring cavity is re- 
alised. It is the following. The microwave power, that is 
microwave photons having occurred in the sample remain 
[are captured and entirely absorbed] in the sample like to 
microwave power in measuring cavity with metal walls. 
Any transmission of microwave power through the samplr 
does not take place. It is the sufficient condition for the 
formation of the second cavity in measuring cavity. In 
other words the sample itself becomes to be the second 
cavity, that clearly was observed to be additional nar- 
row peak on microwave generator generation zone oscil- 
logramme. By AF-mode registration the microwave gen- 
erator frequency follows the changes in the frequency of 
the sample-cavity, determined by time dispersion of the 
sample-cavity. The changes in frequency of the order of 1 
MHz are relatively large and the resonance phonons with 
starting frequency remain in coherent state in the sample 
[they cannot interact immediately with microwave pho- 
tons, which have now the other energy value]. In fact 
lifetime of coherent state of resonance phonons is deter- 
mined by their interaction with nonresonance phonons 
and it can be very long, if given interaction is relatively 
weak. It takes place in the anthracite samples studied 
and it seems to be the consequence of their disordered 
structure. By the sweep back of the static magnetic field 
Hq (^^ < 0) the well known phenomenon of acoustic 
spin resonance (ASR) is realised, which is characterised 
by own Rabi oscillation process, presented in Figure 1. 



It is especially interesting, that acoustic Rabi oscillation 
process has quantum nature with discrete dependence on 
the time (time sweep is connected with the static mag- 
netic field sweep). The possibility to register the ASR 
by means of lOO-kHz modulation canal is determined by 
the following factors. AF unit produces frequency mod- 
ulation of microwave power, which becomes to be not 
srongly monochromatic. At the same time the quantum 
Rabi oscillation process itself is not monochromatic pro- 
cess. Therefore the nonmonochromatic modulation of the 
static magnetic field is taking place, and the Fourier com- 
ponent of given nonmonochromatic modulation is regis- 
tered by 100-kHz modulation canal. Given conclusion is 
confirmed by the polarity of derivative of ASR responce. 
Really, it is substantial, that observed polarity of deriva- 
tive of ESR- ASR responce by usual registration in AF- 
mode does not depend on sweep direction. At the same 
time, when ESR spectra by ^ > and by ^ < 
are tuned for the absorbtion registration the polarities of 
qiven spectra have to be different [1^ . Consequently one 
of the spectrum is emission spectrum. It is evident, tak- 
ing into consideration the energy conservation law, that 
emission, corresponding to ASR phenomenon, is observed 
by < 0. It is understandable, that by the registra- 
tion in CF-mode the process of accumulation of resonance 
phonons, for which the EM-field is the main heat reser- 
voir, will be not realised. It is the consequence of effective 
back process, since resonance phonons have the same en- 
ergy with resonance photons and the only ESR responce 
can be registered in full accordance with experiment, see 
Figure 9. 

Although the emission. Figures 6, 8, is coherent, it is 
not maser effect, since the amplification is absent (maser 
effect seems to be in principle impossible for two level 
system). Then the hysteresis effect (compare Figure 5 
and 7), which is observed by registration in AF-mode, 
can be explained in a natural way like to Stokes lumines- 
cence, where the emission spectrum is shifted to more low 
frequences in comparison with excitation spectrum to be 
consequence of interaction with nonresonance phonons. 
Given conclusion is confirmed by the study of the de- 
pendence of the Stokes shift value in ESR- ASR spectra 
of anthracite sample. Figure 11, and by the study of 
the dependence of the effective linewidth in ESR spec- 
tra for both AF and CF registration modes. Figure 12, 
on the amplitude of 100 kHz high frequency modulation 
of static magnetic field. Really passage velocity through 
the resonance is determined by product of LdmHm, where 
ujm is modulation frequency value. It means, that by 
passage velocity increase the direct interaction of spin- 
photon subsystem with nonresonance phonons [which 
seems to be playing the same role, that the interaction 
with phonons by Stokes shift formation in luminescence 
studies] can be "turned off', then emission and absorb- 
tion spectra have to be registered without Stokes shift. 
Therefore the dependence of the Stokes shift value on 
the amplitude of 100 kHz high frequency modulation of 
static magnetic field has to to have steplike view and 
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Stokes shift has to vanish at high modulation amphtude 
side of given step. It really takes place, see Figure 11. 
Moreover Figure 11 allows to estimate the lifetime of 
coherent state in spin-photon subsystem from the Hm- 
value (in frequency units) in inflection of given curve. It 
gives the value Tg = 1.3 x lO^^s, which gets to typical 
transverse T2-value range for many systems, described 
by classical Bloch equations. The possibility to "turn 
off' the direct interaction of spin-photon subsystem with 
nonresonance phonons by AF-mode registration means 
on the other hand, that the dependences of the effective 
linewidth in ESR spectra by AF-mode registration and 
CF-mode registration have to be qualitatively different. 
In AF-mode at small modulation amplitudes the direct 
magnetic dipole interaction of spin subsystem in coherent 
state between the centers themselves takes place and it 
will lead to broadening of the registered ESR line, at that 
it seems to be evident that linewidth value has to increase 
with the moving off to the left and to decrease (without 
taking into account the modulation broadening) with the 
moving off to the right from the -ffm-value, corresponding 
inflection in Figure 11 value (0.274 G). It will be true, if 
genuine linewidth value is more than 0.274 G. It is evi- 
dent from CF-mode registration, that genuine linewidth 
value is less than 0.274 G, it is equal to 0.18 G. It means, 
that linewidth at inflection (that is at H„i — 0.274G) 
will be modulationary broadened. The competition of 
modulation broadening and magnetic dipole averaging 
processes (like to those ones in liquid systems) shifts the 
minimum position to the value « 0.4G, see Figure 12, 
in AF-mode dependence of the effective linewidth. Tak- 
ing into account the value of modulation broadening at 
« 0.4G, obtained from CF-mode dependence of the ef- 
fective linewidth, we can also evaluate the value of T,. 
It is equal to Tg — 1.36 x 10~^ s. Some difference of 
the value of T^, obtained from experimental dependences, 
presented in Figures 11 and 12, is in limits of experimen- 
tal inaccuracy, estimated to be equal to lO^^s. The value 
of Ts allows to evaluate the range for the time of spin 
polarization transfer, corresponding to so called spin dif- 
fusion time Td by Dyson effect in metals. In the case 
of normal skin effect the value of ^ = 23.4 ± 0.5 means, 
that ^ < 10"'' and consequently Td < 1.3 x 10^^°s. 

The key to structure of ESR centers can be found by 
the analysis of the central part of ESR spectra by regis- 
tration with automatic microwave frequency adjustment, 
presented in Figures 7 and 8. It is seen from Figure 7, 
that Rabi oscillation process is also takes place, although 
spectral resolution is not very good, since the modula- 
tion amplitude is coinciding with linewidth. The most 
interesting is the shape of background line. In both the 
spectra the shape of background line has kink form and 
can be be approximated by tanh-function. It can mean, 
that the model of the ESR centers has to be similar to 
Su-Schrieffer-Heeger (SSH) model of topological solitons 
in trans-polyacetylene [30| . The g-f actor value, which is 
very near to g-factor value of neutral solitons in trans- 
polyacetylene testifies in favour of given proposal. It can 



mean, that structural element of anthracite includes car- 
bon backbone of trans-polyacetylene. At the same time 
anthracite is structurally disordered material. It means, 
that the solitons has to be pinned and they have to give 
usual ESR spectra. To explain the appearance of kinks, 
at that it is substantionally, that they are characterized 
by different parameters for absorption and emission spec- 
tra, and the possibility to register Rabi oscillations, we 
have to consider joint {EM-field 4- magnetic dipole -I- 
resonance phonon field} system. EM-field and resonance 
phonon field produce radiation communication in the di- 
rection of EM-field propagation and in time. It means 
that we have ID chain, at that the chains of two kinds 
are produced in time - with the sequence - absorption 
time segment - emission time segment and with reversal 
sequence. It is substantionally, that given time segments 
are different in their values. In fact the dimerisation take 
place like to trans-polyacetylene, however instead space 
z-axis along t-axis. It is in agreement with known equal- 
ity of rights of coordinates in Minkowski space. In other 
words the time has to be considered being to be quantized 
in given case. The junction of given time ordered chains 
produces time topological violations of two kinds. The 
first topological violation can effectively absorb photons 
and emit resonance phonons, the second topological vio- 
lation has reversal properties. So, we come to the model 
of instanton, which is like to some extent to SSH-soliton 
model. It is understandable, that the instantons of both 
kinds can free, that is without additional activation en- 
ergy, move along t-axis like to SSH-solitons along space 
z-axis. The very short value of Td < 1.3 x lO'^^s is 
strong experimental support in the favour of the model 
proposed. Given model allows to explain the rather nar- 
row ESR lines in structurally disordered materials like to 
sample studied, at that they are in more than 50 times 
narrower in comparison with ESR lines in very good or- 
dered related material - carbon nanotubes, produce d by 
hig h energy implantation of diamond single crystals jll | , 
|12| . The instantons seems to be two level noninteracting 
between themselves objects. Consequently they charac- 
terised by discrete energetic structure and the splitting 
into energetic band does not taking place, so we see that 
instanton model explain the appearance of narrow lines 
in anthracites in a natural way. It is the first experi- 
mental detection of instantons (to our knowledge) . The 
known instantons in the literature are considered to be a 
special kind of vacuum oscillations in gluon field theory 
and also in gravitation field theory [3]| . They have quite 
another origin. 

Therefore we have observed new kind of coherent 
states. They are modes of resonance phonon field. The 
coherence time value in the samples studied is evaluted 
to be equal at room temperature ~ 10ms, at that at rel- 
ativelly low microwave power and it will increase with 
power increase, achieving the maximal known values at 
present and even more. Substantial advantage for pos- 
sible applications is that, that the responce on external 
EM-field is very strong. At the same time the detection 



12 



of single spin is rather technically difficult task. The H- 
field dependence allows to separate the phonon modes, 
the dependence on H-field direction seems to be also es- 
sential for practical applications. 

III. CONCLUSIONS 

QED - QPFT model for dynamics of spectroscopic 
transitions in ID multiqubit exchange coupled system 
is g eneralized by taking into account the earlier proof 
[14{ , that spin vector is quaternion vector of the state of 
any quantum systen in Hilbert space defined over quater- 
nion ring and consequently all the spin components has 
to be taken into account. New quantum phenomenon 
is predicted. The prediction results from the structure 
of the equations derived and it consists in the following. 
The coherent system of the resonance phonons, that is, 
the phonons with the energy, egualed to resonance pho- 
ton energy can be formed by resonance, that can lead to 
appearance along with Rabi oscillations determined by 
spin (electron)-photon coupling with the frequency fl^^ 
of Rabi oscillations determined by spin (electron)-phonon 
coupling with the frequency Vl^^^ . In other words QED 
- QPFT model predicts the oscillation character of quan- 
tum relaxation, that is quite different character in com- 
parison with phenomenological and semiclassical Bloch 
models. Moreover if absolute value of electron-phonon 
coupling constant | | , which characterisires the interac- 



tion with S'^- and components of the spin of jth chain 
unit is less, than electron(spin)-photon coupling constant 
5, the model predicts that the second quantum Rabi os- 
cillation process will be observed by stationary state of 
two subsystems {EM-fied -I- magnetic (electric) dipoles}, 
and it will be registered in quadrature with the first quan- 
tum Rabi oscillation process. The second quantum Rabi 
oscillation process is governed by the formation of the co- 
herent system of the resonance phonons. Therefore along 
with absorption process of EM-field energy the coherent 
emission process takes place. Both the quantum Rabi 
oscillation processes can be time-shared. For the case 
of radiospectroscopy it corresponds to the possibility of 
the simultaneous observation along with (para)magnetic 
spin resonance the acoustic spin resonance. The second 
(acoustic) quantum Rabi oscillation process can be de- 
tected even by stationary spectroscopy methods. The 
model of new kind instantons is proposed. All the predic- 
tions are experimentally confirmed for radiospectroscopy. 
It has been found in particular, that the lifetime of coher- 
ent state of collective subsystem of resonance phonons in 
anthracites of medium-scale metamorphism is very long 
and even by room temperature it is evaluated in ^ 10 
ms. The phenomenon of the formation of the coherent 
system of the resonance phonons can find the number 
of practical applications, in particular it can be used by 
elaboration of logic quantum systems including quantum 
computers and quantum communication systems. 
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